function gamma=computeGamma(identityC, powersAB)
gamma = identityC * powersAB;
end